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Abstract 

We discuss QCD evolution equations for two and three particle correla- 
tion functions of quarks and gluon fields in a hadron which describe devel- 
opment of the momentum distribution of a parton system with a change of 
the wave length of a probe which resolves it. We show in a general case of 
two-particle correlators how the four-dimensional conformal algebra and the 
known pattern of conformal symmetry breaking in QCD can be used to solve 
the complicated mixing problem of local operators under renormalization and 
compute economically anomalous dimensions of quark and gluon composite 
operators. An extension of QCD to Af — 1 super Yang-Mills theory and use 
of superconformal anomalies arising after quantization allows to derive non- 
trivial relations between the anomalous dimensions. For three-parton systems 
the conformal symmetry alone is not enough to solve the three-particle prob- 
lem. We show that in milticolor limit of QCD there arises an extra conserved 
charge describing the solitonic motion of the system of particles. The problem 
admits a one-to-one correspondence with certain spin chain models which are 
exactly solvable. 
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Wc discuss QCD evolution equations for two and three particle correlation func- 
tions of quarks and gluon fields in a hadron which describe development of the 
momentum distribution of a parton system with a change of the wave length of a 
probe which resolves it. We demonstrate how broken space-time and hidden QCD 
symmetries serve to understand the structure and ultimately solve these equations. 

1 Unraveling layers of matter: from atoms to partons 

A wisdom goes back to ancient Greeks who philosophized that the matter 
consists of tiny particles — atoms. However, the atomic strticture remained 
a puzzle till the beginning of the 20th century when the radioactivity was 
discovered and used by Rutherford in his seminal experiments on large angle 
scattering of a particles off atoms which suggested that the atom bears a lo- 
calized core — the nucleus. On the other hand, electron beams were found 
to pass through atoms with no or very little deflection forcing Lenard to hy- 
pothesize that atoms have wide empty spaces inside, a particles having much 
larger wave length comparable to the nucleus size scattered more frequently 
with a low intensity beams than their 'cousins' /3 particles having much smaller 
wave length and which, for available luminosity, had a very low cross section. 
Similar experiments but rather with light sources or neutrons are exploited 
nowadays to study the crystal structure. If one puts a crystal in front of a 
source of visible light, the object just leaves a shadow on a screen behind it 
and one does not see elementary building blocks which form it. Obviously, the 
visible light, having the wave length A^ght ~ 0.4 — 0.7 /itm cannot do the job 
and resolve the internal structure of a crystal. The size of an individual atom, 
say hydrogen, is of order Tatom ~ (acmWe)^^ ~ (10 KeV)""^. Therefore, to 
'see' atoms in crystals one has to have photons with a comparable or smaller 
wave length < ratom, or equivalently, of energy > r~^^^. To do this 
kind of 'nano-photography' one need a beam of X-rays. To go beyond and 
look into the structure of atoms one needs even smaller wave lengths. 

After the discovery of nucleus' building blocks — nucleons, i.e., protons 
and neutrons, the attention has been shifted to the study of these 'elemen- 
tary' particles. However, their elementarity has been questioned starting from 
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Figure 1. Partonic picture of inclusive (a) and exclusive (b) scattering of a 7*-probe from 
the nucleon. 



Stern's measurements in 1933 that demonstrated a large proton's magnetic 
moment. Hofstadter's experiments with elastic electron scattering off nucle- 
ons, eiV — > eW, and measurement of their charge distribution revealed that 
the nucleon has a spatial extent of order 1 fm. To probe femtometer scales 
one has resort to scattering experiments with high energy beams. Inelastic 
lepton- nucleon scattering experiments, eiV — > e'X, conducted at SLAC ful- 
filled this goal and led finally to the discovery of a new layer of matter by 
observing events with the transfer of a large momentum from the electron to 
the proton. If the latter would be a hard ball or a diffuse distribution of mat- 
ter, such kind of scattering would be improbable. Rather it was explained by 
conjecturing the existence of point-like constituents inside the nucleon which 
absorb a highly virtual 7*-quantum emitted by leptons, see Fig. [I] (a). Anal- 
ogous representations holds for the elastic scattering alluded to above. Fig. |l] 
(b). These partons are nothing else than QCD quark and gluons described by 
a rigid SU (3) gauge field theory. 



2 Hard scattering and QCD 

A QCD scattering amplitude is given by a successive sequence of interactions 
and free propagations through rife intermediate states created by QCD field 
operators. At rtth order in the coupling g = ^/Amxl it reads in the time- 
ordered form 

J dtn j dtn-1 ■ ■ ■ j dti Gnjn-l{Ui)Vn-lGn-l,n-2{tn-l) ■ ■ -VlGlfiiti) . 

— oo — oo — oo 

The summation in repeated indices runs over states which are eigenfunctions 
of free equations of motion. Vn^m — {n\V\m) is a matrix element of a pertur- 
bation, — the interaction potential, — and the free particle propagation is 
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expressed in terms of a time-underintegrated free action, S, which defines a 
Hamiltonian of the noninteracting system, 



At this point, one immediately concludes that there is no sensitivity to long- 
time acales as long as the total phase, the free particle action is time depen- 
dent In this case, cancellations of long-time contributions occur due to 
oscillatory nature of the integrand. However, if the phase is stationary, which 
corresponds to a stationary classical action, particles travel along their classi- 
cal trajectories and the amplitude becomes sensitive to large-time dynamics. 
Since at large space-time scales as is large, one cannot evaluate the amplitude 
in a perturbative expansion and effects of quark confinement are relevant. In 
the opposite situation, which goes under the name of infrared safe, an ob- 
servable can be computed as a power series in QCD coupling alone. These 
include the total inclusive atot and jet crjct cross sections in electron-positron 
annihilation, etc. Indeed in atot ~ {^\T{jj,,jfj.}\0) , once a pair of quark and 
antiquark is created by the source , they travel back-to-back with a speed of 
light and cannot reassemble back into a physical state absorbed by jj^. Since 
there are no classical trajectories which allow this processes, the observable is 
not sensitive to infrared physics according to the criterion spelled out above. 

A hard scattering cross section having at least one hadron in the initial 
state cannot be infrared safe since in a preparation of the asymptotic hadron 
state partons which form it had interacted strongly with each other for a long 
time in a bound state. Obviously, the wave functions of quarks in a bound 
state differ from what one would have if they were free. So this inevitably 
affects the cross section. Let us discuss in detail the deeply inelastic scat- 
tering example in Fig. ^ (a). The incoming lepton fluctuates into a lepton 
and a photon, e — s- e'7*, and the latter interacts with partons in the target 
|A^) fragmenting into a number of hadrons in the final state \n). Since the 
electromagnetism is very weak as compared to the strong interaction, one can 
restrict the analysis to a single-photon exchange with spacelike momentum 
q, = —Q^. In order to initiate a hard scattering, the electron has to pass 
close to one of the partons, i.e., at the distance ~ to exchange a 

photon of virtuality Q^. This process is described by the amplitude (n|j^|A'^) 
with the local quark electromagnetic current j^. The measurement is totally 
inclusive with respect to final states and only the scattered lepton is detected. 
The cross section reads 



Gn^m{t) = exp{it {Sn - S,n)} ■ 




(1) 
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Let us turn to the physical picture of the deep-inelastic event in the photon- 
hadron center-of-mass frame. For a highly virtual photon, the transverse 
distance probed by it, in a Lorentz contracted hadron, is of order bz\^ ~ 
Due to a Lorentz dilation, the virtual photon 'sees' the partons in a frozen 
state during the time of transiting the target. Therefore, since the number of 
parton which carry the bulk of the hadron momentum is small, the photon 
will 'see' mostly only one parton per collision. The probability for multiparton 
correlations is suppressed by the hard momentum, 

n-parton probability - ((<5z_l)^ /^-l)" 1/ (Q^^i)" , (2) 

where = Tri?^ is the transverse area of the nucleon of the radius i?jv. 
These power suppressed multiparticle correlations go under the name of higher 
twists. Thus, to leading approximation one can restrict considerations to a 
single-parton scattering at high . As compared to the inclusive annihilation 
mentioned above, the underlying physical picture for the forward Compton 
scattering admits a classical trajectory. A quark taken from the hadron ab- 
sorbs the virtual photon and, as a result, accelerates. Then it reemits a photon 
and falls in the same momentum state. After the energy is freed into the final 
state the parton merges back into the parent hadron. As we already discussed 
above, the process is not infrared safe and depends on the quark binding in- 
side the nucleon. The higher is the virtuality of the virtual photon Q^, the 
shorter are the distances traveled by the parton. The points of absorption 
and emission are separated by a light-like distance. The character of relevant 
distances in the Compton amplitude is a consequence of deep Euclidean kine- 
matics oo. Large and energies p = p ■ q, at fixed a;B — Q^/(2i/), 
probe short distance and time structures of the process, respectively. To iden- 
tify scales involved in deeply inelastic scattering let us go to a reference frame 
where the target proton is at rest and the virtual photon's three-momentum 
is directed along the z-axis, q — {q'^, 0, 0, q'^). Then 

For light-cone components of the momentum transfer, q± = (go i (73) /a/2, we 
have for a large momentum transfer 

g_ - QV (Mxb) , q+ ^ MxB . (4) 

The integrand in Eq. (|l]) is an oscillatory function and thus gives vanishing 
result unless the distances involved are 

z_ - 1/ (Mxb) , z+ - Mxq/Q^ . (5) 
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Figure 2. QCD amplitude for deeply inelastic scattering from a quark in the nucleon and 
resulting scaling violation (evolution) of the probability to find a parton the with a given 
momentum fraction x of the parent hadron with a change of the virtuality of 7* . 



Therefore, provided transverse separations z± are small, the deeply inelastic 
scattering probes strong interaction dynamics close to the light-cone ~ 0, 
and we can neglect the dependence on all coordinate components except for 
Z-. The latter is called loffe time and has the meaning of the longitudinal 
distance probed in the process. Its Fourier conjugated variable is the momen- 
tum fraction x of the nucleon carried by a parton interacting with the probe. 
In lowest order approximation x — x^. 

A struck quark can emit a gluon or many gluons before actual interaction 
with the photon, see the left drawing in Fig. |[ When it spills a gluon the 
scattered quark looses a portion of its original momentum. Although it can 
be small for a single event, quarks tend to emit more being a relativistically 
moving particle. Thus, the bremsstrahlung of many gluons drifts the quark 
into the low momentum fraction region. Since gluons in turn can decay into a 
quark-antiquark pair, there is a proliferation of partons with small momentum 
fractions and correspondingly a decrease in the large- a; domain. Thus, the 
probability to find a quark with small momentum fraction is higher at larger 
Q^, see the right drawing in Fig. |^. The same phenomenon can be viewed 
from the resolution perspective discussed above. The higher the smaller 
is the distance probed by the photon. Therefore, one sees more and more 
partons in the cloud forming the parent 'fat' quark. 

Since the hard subprocess occupies a very small space-time volume, the 
scales involved in the formation of their nonperturbative wave function are 
much larger, of order of a typical hadronic scale, 1 GeV. Hence, it is quite 
likely that they are uncorrelated and will not interfere. Thus, although the 
process depends on the hadronic state the parton has come from, this is 
basically irrelevant for the hard interactions occurred. Moreover, all final state 
interactions cancel in the deeply inelastic process. This is exhibited by the 
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relation (|l|) stemming form the optical theorem. Thus, there is no sensitivity 
to soft final-state interactions. The quantum mechanical incoherence property 
of physics at different scales results into the factorization property of the cross 
section (|l|) . This is the base for the predictive power of perturbative QCD in 
many hadronic processes which one writes as 

a {xB,Q^) ^ J dxC (xb/x; Q^fi^) F (x; ^i^) , (6) 

where is a parton distribution, — the probability to find a parton in the nu- 
cleon with the momentum fraction x, — and C is a perturbatively computable 
short-distance coefficient function. 



3 Evolution equations 

Due to intrinsic ultraviolet divergences, parton distributions and coefficients 
functions in Eq. (||) depend on a renormalization scale However, it can- 
cels, as it has to, in the total result since the left-hand side being a physical 
observable does not depend on the arbitrary scale /Lt^. Thus, factorization 
implies evolution as a consequence of this independence, 

0. (7) 



d In /i^ 

Separation of variables results into complementarity equations for F and C 

■^^Cix; Q^fi^) = - j dx' P{x/x'-as)C{x'- QV/x^) , 

J dx' Pix/x';as)F{x'; ^i^) . (8) 

The coefficient function C is computable perturbatively and so is the evolution 
kernel P. The second of these equation is an inclusive evolution equation and 
allows to predict _F at a scale from the function F determined at some 
reference scale fiQ. Obviously, a change in fi^ is in one-to-one correspondence 
to a change in ^ 

In the previous discussion we talked only about quark momentum proba- 
bilities. To describe a real experiment one has to account for another type of 
partons which are equally important dynamical degrees of freedom — gluons. 
Gluons contribute on equal footing with quarks and both enter in a multiplet, 



"^We will not discuss presently the evolution in longitudinal momentum with fixed transj«rse 
scales. This domain is described by another type of evolution equations addressed in 0. 
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Figure 3. Elementary a — » fee parton splitting amplitudes encountered at each step of the 
evolution (left) and a fit to experimental data □ (right) implementing them (dressed with 
next-to-leading order perturbative corrections) in the evolution equation (p|). 



(FqjFg). This multiplet obeys a matrix evolution equation of the type (| 
with a quark-giuon mixing matrix of kernels 



qqp lap 
gqp gap 



with -'Fix) ^ \Aa^bc\\ (9) 



where Aa-*bc represent amplitudes of elementary splittings, demonstrated in 
Fig. p. The solution to the evolution equation (||), where one includes also 
QCD radiative corrections as well, results to a perfect description of experi- 
mental data for the unpolarized deeply inelastic cross section, shown in the 
same figure. 

We consider below a more general case of functions which interpolates 
between inclusive parton distribution and exclusive distributions amplitudes 
and, thus, possess a very rich structure. They are known as generalized parton 
distributiojis. F{x,ri) and are measurable in a number of hard exclusive 
reactions Qfl'ElllS such as in Fig. ^ (b) with a photon emission into the final 
state. They admit an operator representation which reads, e.g., for quarks 

{N{p2M-z^)q{z^)\N{pi)) = f dxe-'"'-^P'+P^'>+F{x,Tj). (10) 
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It is a correlation function of fields on the light cone separated by the distance 
2z_ which acquire a scale dependence under renormalization. We omitted an 
inessential to our present consideration Dirac matrix projecting quark spinor 
indices. Due to different momenta of incoming and outgoing hadrons, the 
specifics of these functions lies in a presence of a non-zero ^-channel momentum 
transfer rj = {p2 —pi)+/2, varying which one probes correlations of partons in 
different momentum states. The evolution equation has the same schematic 
form as (||), 



F{x,-q;^?)^ / dx' K{x/tj,x' /T];a,)F{x',7j;fi'^) . 



din /i^ 

The evolution kernel interpolates between two known limits, i.e. 



(11) 



lim K{x/ri,x' /ri;as) = P{x/x';as) , lim K{x/ri,x' /r]]as) = V{x,x']as) , 

the fully inclusive case we have discussed above (77 — 0) and another one 
which corresponds to distributiaLi,,amplitudes which enter the description of 
a number of exclusive processes E3I13. 

The moments of the function F{x, rf) correspond to the operator matrix 
elements of local operators, e.g., for quarks. 



dxx'Fix,!]) = {P2\qd'+q\Pi) , 



(12) 



where d = d — d ■ Under renormalization the operators in Eq. ( |12D get mixed 

with operators containing total derivatives d = d + d but, indeed, the same 
total number of derivatives as a consequence of Poincare invariance 



n,k=qdld'+q. 



(13) 



Since the operators with total derivatives matter for each moment we have to 
diagonalize a j x j matrix equation. 



(14) 



with zeros below the diagonal. The matrix has obviously a triangular form 
since the operator containing at least one total derivative cannot mix with 
operators involving none of them. Otherwise one would run in a contra- 
diction for the forward matrix elements where total derivatives are irrelevant 
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{p\d{- ■ ■)\p) — 0. Another important observation is that the diagonal elements 
of the anomalous dimensions are the moments of the splitting functions 



ab D 

h 



-2 J dxx^^^Pix). (15) 



4 Use of conformal symmetry 

Recall the diagonalization of the Schrodinger equation for a hydrogen atom. 
The 0(3) symmetry of the potential allows to reduce the three-dimensional 
problem to a one-dimensional one, i.e., one performs a partial wave decom- 
position in terms of spherical harmonics. The spherical harmonics are eigen- 
functions of the quadratic Casimir operator of the rotation group. 

For the case at hand one can use the space-time symmetries of QCD 
Lagrangian to simplify significantly and ultimately solve the problem. Classi- 
cal QCD is Poincare, i.e., Lorentz AA^i, and translation V^, invariant as any 
realistic field theory describing nature, a well as due to absence of any di- 
mensionful parameters for massless quarks, invariant under dilatation T) and 
special conformal /C^ transformations tj. Coordinates transform as follows 
under them, 

M ^1 ^ \ 2 _ ^ 

The special conformal transformation can be visualized as a sequence /C^ = 
IV^I of inversion Tz^ = z^j z'^, translation and inversion. This fifteen trans- 
formations form a group, 5*0(4,2). It is easy to construct tha.representation 
of the group on field operators using induced representations EJ. For instance, 
S^q = i[q,Vfj]- — —dfj^q, etc. The eigenvectors of the quadratic Casimir 
operator of the projected coUinear conformal group, 

J2 = ^V+IC^ - i (P + M-+ + 2i) {V + , (16) 

in the basis of the bilinear operators (|l^) are called conformal operators with 
eigenvalue^ 

jl^Oji = JabiJab - l)Ojl , Jab = J + Ja + Jb ■ (17) 

For quarks j^==_L and conformal operators are given in terms of Gegenabauer 
polynomials B'EI (/ > j) 

0,i^qd\c]'''"^ (d+ld + )q, (18) 

''The indices attached to stand for the conformal spin ja = {da + Sa)/2 of fields, having 
the dimension da and spin Sa, which build Oji- 
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which is a very specific combination of local operators TZ jk (|l^). 

Generally, the renormalization group equation reads for the conformal 
operators 

'^i=~^i:iifs,,+j]i^)o,,. (19) 



^ fe=0 

As is well known a renormalization group equation is a Callan-Symanzik equa- 
tion, i.e., the scale Ward identity for a Green function with a (conformal, for 
the case at hand) operator insertion, 

g{zi,...,ZN) = I DqDqDAe'I''''^^'^Ojiq{z,)...q{zk)...A{zN). (20) 

The anomalous dimensions of operators arise due to development of an 
anomaly in the trace of the energy- momentum tensor, 

S^C = -0,, = ^{GI^)\ (21) 

where Pd is a d-dimensional QCD beta function f3d = f3 + g{d — 4)/2. We 
did not display other terms which could not affect the physical sector Ii3. The 
anomalous dimensions show up in the renormalization of an operator product, 

0,i 1 1 d'^z 0^^ [z) = {I'^Sjk + iT) Ou ■ (22) 

Since the conformal algebra is fulfilled on interacting (Heisenberg) field 
operators one can derive a constraint on the anomalous dimension from it. It 
is provided by the commutator 

[P,/C^] =z/C„, (23) 
which results into a commutator equation 00 



9 



= 0, (24) 



in terms of ag-independent matrices a and 6; QCD beta-function P and a 
special conformal anomaly 7^ which arises, similarly to Eq. (^), in the renor- 
malization of the product 

0,1 j d'^z S^C{z) = ^2 l'ikOk,i-i . (25) 

fc=0 
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The special conformal variation of the action is given in terms of a coordinate 
moment of the energy-momentum tensor, 5^C{z) = — 2z_6'^^(z). 

One can immediately draw important conclusions from Eq. (p^). If one 
is interested in the one- loop anomalous dimension matrix, 0{as), one has to 
keep the Ofg-independent term a only and write [7°(°) -|- ^nd(o)^ _ Since 
the matrix a is diagonal 113 we conclude that ~ Q ■ as- At next-to- 

leading order, we have to account for (3 and 7"^ in one-loop approximation. 
This provides via Eq. ( p4| ) an expression for two-loop nondiagonal elements 



7 



ND(1) 



in terms of 7 



D(0 



7 



ND(1) ^ 



/3o and 7'=("\ 
h 



D(0) 



a V 



vD(0) 



c(0) 



(26) 



No two-loop calculation has to be done! The diagonal anomalous dimensions 
7^^ coincide with the anomalous dimension of local operators without total 
derivatives and are available in the literature on deeply inelastic scattering. 



5 Use of superconformal symmetry 

As we shall see momentarily, supersymmetry provides relations between the 
anomalous dimensions of conformal operators, while the extension to super- 
conformal case results into identities between the special conformal anomalies 
7*^ of quark and gluon operators. 

If one puts the QCD quarks into the adjoint representation of the color 
group and sets the number of flavors to one, one recovers M = \ super- Yang- 
Mills theory in the Wess-Zumino gauge. In four-dimensional space-time the 
Lagrangian is invariant under superconformal group, i.e., S'0(4, 2]_iiiacussed 
above plus translational Q and superconformal S transformations EjO'lS, 

-5^9" = ^G^,.^M^e-*^"75e , 5^Bl = -zCi^q^ , 5^D- = 7^759'^ , (27) 

with non-dynamical auxiliary field D°' and transformation parameter ^ = ^0 ~ 
i ;6S,i. The parameter ^0 (Ci) parametrizes T — Q (T — S) transformations. 

One can construct an Af — I chiral superficld from certain linear combi- 
nations oLbpaonic, i.e., qq and gg, and fermionic, i.e., gq, conformal operators, 
see Refs. oEB. So that under ( p7| ) they transform as Wess-Zumino multiplet. 
This covariance simplifies the derivation of relations between the anomalies. 

We follow the same procedure as in the previous section by using the 
commutator algebra of the supeconformal group for the Green function 
G{zi, . . . , zjv) (^0|) and use Ward identities to find out consequences. Because 
of the quantization which require a regularization procedure, the variation of 
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the action develops an anomaly in d-dimensions for the 5-transformation but 
not for the Q-transformation 

6QC = 0, S^C = ±^^,^a,,G;,q^ . (28) 

Although there arise contributions in both due to explicit supersynimetry 
breaking by gauge fixing, however, these sources cannot affect gauge-invariant 
quantities we deal with. For our purposes we need two commutators from the 
algebra 

[Q,2?]- = ^Q, [S,/C^]_-7m5. (29) 

In view of Eq. ( p8| ) the first one is non-anomalous and results into a 
PjUabejr of relations between the elements of the quark-gluon mixing matrix 
ll3c3'Eil, with the most trivial example involving only the diagonal elements 

ii^f + J '"if = I ""if + ""if ■ (30) 



The second commutator in Eq. (|29|) possesses the superconformal anomaly 
( p8| ) and, as a result, relations of the type (^, but with 7 being replaced by 
7^^, are affected in a controllable manner by a computable addendum A which 
arises from renormalization of the operator product E3 

0,1 i jd^z 5^C{z) - ^2kOk^-i ■ (31) 

fc=0 

The knowledge of a complete set of such relations allows to compute all 
entries of the quark-gluon mixing matrix from the knowledge of only one 
element, e.g., '?'?7° -|- 997^0^ ^j^jg procedure has successfulliplpeen apphed in 
the reconstruction of two-loop off-forward evolution kernels c^. 

6 Integrability of three-particle problem 

As we explained in the introduction, multi-particle correlations arise as power 
suppressed effects in hard cross sections and exhibit a genuine quantum me- 
chanical interference of hadron wave functions with different number of con- 
stituents. The first non-trivial and the most interesting examplej-is three- 
particle functions which emerge in hard forward Ej and off-forward l3 Comp- 
ton scattering in the form of quark-gluon-quark {qgq) and three-gluon {ggg) 

arrelations. Another prominent example is baryon distribution amplitudes 
. For the sake of definiteness, let us discuss quark-gluon-quark correlators. 
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Figure 4. Large- A^c reduction of three-particle qgq Hamiltonian. 



As all other twist-three sectors, they are distinguished from all even higher- 
particle sectors by the fact that their basis can be reduced to the one consisting 
only of quasi-partonic operators E3. This implies that the three-particle kernel 
of the leading order evolution equation in QCD coupling constant, 

-^F (W; ^,^) = ^lf[ dx, 5 X, -^K {{x}- {x'}) F {{x'}- f,') , 

(32) 

with {x} = xi,X2, X3, is merely given by the sum of pair-wise interactions of 
twist-two kernels discussed in the previous sections, K = Kqg-\-Kqg-\-Kqq. As 
a result of the conformal covariance, established there, the pair-wise kernels in 
the basis of local conformal operators, to be called Tiab later, depend only on 
two-particle Casimir operators, J^^ = JabiJab — 1)7 of the conformal group, 

T-C = Tiqg (Jqg) + Hqg (Jqg) + Hqq (Jgq) . (33) 

The color structure of pair- wise interaction allows to simplify further the diag- 
onalization problem by noticing that Tiqq vanishes as l/iV^ relative to 'Hi^q^q-^g 
in multicolor limit, i.e., Nc —* 00. Therefore, the three-particle problem on 
a ring is reduced to the one on a line with the interaction of end sites being 
neglected, see Fig. ^. Our cherished goal is to solve a Schrodinger equation 

H'^'F EE {n^g + n^^g) F = NcS'^F . (34) 

The explicit form of the pair-wise Hamiltonians depends on the i-channel 
quantum numbers. For instance, for the operators of the Dirac-color structure 
qapif.fGljpq they are 

^^abiJab) = i> [Jab + 3/2) + i, [J^b " 3/2) - 2^{l) , (35) 

where ^{x) = din F (x) / dx and the indices run over a — {q,q) and b = g. The 
large-iVc Hamiltonian 7i°° describes a three-site inhomogeneous integrable 
open spin chain model. Integrability means that it possesses a complete set of 
integrals of motion matching the number of degrees of freedom. The problem 



evolution: submitted to World Scientific on February 1. 2008 13 



admits a representation in terms of the 7?,-matrix acting on a_tensor product 
of vector spaces Va ®Vb and satisfying Yang-Baxter relation l3, 





One of the solutions to the Yang-Baxter equatiDn is given by TZabW 
X)/r{Jab — The construction by SklyanintZl allows to find the H; 
corresponding to this 7^-matrix via the following graphic equation 

2 dX 



(36) 








q S q 




ln< 


'1 


g g 
g g 


II 


A=0 









(37) 



The spectral parameter A in the q and q sites is shifted as A ± 3/2. The 
picture is self explanatory: the three 7?,-matrices with dimensions of spaces as 
displayed are aligned next to each other, then reflected with the unity matrix 
and traced in the auxiliary space with their inverse. The resulting object 
is called the transfer matrix. Its logarithmic derivative determines the qgq 
Hamiltonian up to an additive constant. 

The Yang-Baxter equation holds for an arbitrary dimension of the vector 
spaces. When the auxiliary space is two-dimensional, one finds another solu- 
tion to (p^), 'R-a,i/2W = Al -I- cr* J*. Calculating the transfeK-mAtdx for the 
two-dimensional auxiliary space, one identifies two operators L^tSnJ 

J' - Jl + + A, ' 2 = {^9^ Ag\ + - 2(3/2)^ {J\, + Jl) , (38) 

which commute with the Hamiltonian (^7|) by construction since all of them 
are deduced from bundles satisfying the Yang-Baxter equation. 

Thus instead of diagonalizing the Hamiltonian one can attempt to solve a 
much simple eigenvalue problem for the fourth-order differential operator Q. 
Because of commutativity, the Hamiltonian shares the same eigenfunctions 
Thus the eigenvalues of 7i°° arc functions of quantized values J and Q of 
the integrals of motion and Q. The problem was solvecLusing WKB, i.e., 
large- J, approximation with the result for energy levels P^rr^ri 

£^{J) = V'(J + 3) + ?A(J -I- 4) - 2?A(1) - 1/2 , (39) 



\J,Q) = 2 In J -I- 25fte (3/2 -I- i^{2Q/J'^ - 3) /4) - 4i/^(l) - 3/2 , 



where the special value is separated from the rest of the spectrum by a 
gap and was known for some time e2I. Full ap:pupt of this and other exactly 
solvable cases can be found in the literature E3E3c3. 

We would like to thank the organizers of the workshop for such an inter- 
esting meeting and INT (Seattle) for support. This work was supported by 
the US Department of Energy under contract DE-FG02-93ER40762. 
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